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V~j ' Abstract 
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Rare radiative decay — rjTT^'y was studied with SND detector at VEPP-2M 
electron-positron collider and its branching ratio was measured: B{(f) rjir^'y) = 
(0.88 ± 0.14 ± 0.09) ■ 10"^. Significant contribution of the ao(980)7 interme- 
diate state was observed in the decay. The result is based on total integrated 
luminosity corresponding to 2 ■ 10^ produced (j) mesons. 
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Introduction. The first obser- 
vation of the rare radiative decay 

r77r°7, (1) 

and measurement of its branching ratio 
were performed in Novosibirsk by SND 
detector at VEPP-2M e+e" coUider 
0. The analysis was based on data col- 
lected by SND in 1996 H- Later the re- 
sults were confirmed by CMD-2 group 
in their recent publication Q. Results 
of the present work are based on analy- 
sis of all SND data collected in the vicin- 
ity of meson in 1996-1998 0. 

Reaction is especially interesting 
in connection with the scalar ao(980) 
meson problem, which is being actively 
discussed in the literature. At present 
there is no generally accepted viewpoint 
on the nature of Oq, its quark struc- 
ture is still not well established, and 
several models exist including modifi- 
cation of qq scheme 0, KK molecu- 
lar model 0, and 4-quark model P). 
It was suggested in that the decay 
(f) — > ao (980)7 ~^ rjTT^'-f may serve as 
a probe of the ao-meson quark struc- 
ture. Theoretical predictions for the de- 
cay branching ratio vary from 10~^ for 



the simple two-quark and KK molecu- 
lar models up to 10~^ in the 4-quark 
model [§, |l^. There also exist some 
models in which high rate of the decay 
(IID can be achieved without assumption 
of 4-quark structure of the ao meson 
[ pj] ] . Another possible mechanism of the 
(j) ?77r°7 decay is (p p^ir^, p — *■ 
?77. The vector meson dominance model 
prediction for this branching ratio is 5 • 
10-^ [g g. Detailed study of the ^ 
rjn^'-f decay may provide decisive infor- 
mation on ao(980) meson problem. 

Experiment. The SND m is 
a universal nonmagnetic detector. Its 
main part is a 3-layer electromagnetic 
calorimeter consisted of 1630 Nal(Tl) 
crystals. The energy resolution of the 
calorimeter for photons can be described 
as aE/E = A.2%/^E{GeV) 0, the 
angular resolution is close to 1.5°. The 
solid angle coverage is 90% of An stera- 
dian. 

The data used for the study of — *• 
rjn^'y decay were collected in 1996-1998 
0. Nine successive scans of the energy 
range 980-1040 MeV were performed. 
The data were collected at 16 beam en- 
ergy points. The total integrated lumi- 
nosity in the experiment is 12 pb^^ and 



total number of produced mesons — 
2-107. 

Event Selection. Main sources 
of background for the process under study 

e~^e~ — > — > ?77r'^7 — > 57 (2) 

are the following 0-meson decays: 

e'^e~ — > — > n^n^'j 57 (3) 

e"*"e~ ^ ^ ?77 ^ 37r°7 7j, (4) 
e'^e" ^ ^ KsKl — > neutrals (5) 
and a nonresonant process 

e^e^ — s> CJTT*^ ^ 7r°7r'^7 — > 57. (6) 

The process (H) does not produce 67 
events directly but can mimic the pro- 
cess @ due to either merging of close 
photons or loss of soft photons through 
openings in the calorimeter. The pro- 
cess (U) contributes due to Ks — ^ 7r°7r° 
decay accompanied by either nuclear in- 
teraction of the Kl meson or its decay 
in flight. 

Primary event selection was based 
on simple criteria: the number of recon- 
structed photons is equal to five, there 
are no tracks in the drift chamber, the 



total energy deposition Etot ranges from 
0.8 up to 1.1 of the center of mass en- 
ergy 2Eo, the total transverse momen- 
tum of photons is less than 0.15Etot/c. 
In order to suppress background from 
the processes (H) and (|]) a special "pho- 
ton quality" parameter ( [0] was used. 
For i-th reconstructed photon the Q is 
a minus logarithm of likelihood for the 
corresponding transverse energy depo- 
sition profile observed in the calorime- 
ter to be produced by a single photon. 
For multiphoton events ( is defined as 
a maximum Q. The cut ( < sup- 
presses the background from the process 
by a factor of two, reducing the de- 
tection efficiency for actual 5-7 events 
by only 8%. To suppress beam back- 
ground photons which appear mostly in 
the calorimeter areas closest to the beam 
and are relatively soft, additional cut 
was imposed on polar angles of the two 
softest photons in an event: 32° < 6^4, 6^5 < 
148°. 

For events which passed the cuts de- 
scribed above kinematic fitting under two 
alternative hypotheses was performed and 
corresponding values of calculated: 

• an event is one of the process e^e^ — 



57; the value is denoted as Xs^; 

• an event is one of the process e^e' 
37 with two additional stray pho- 
tons; the value is denoted as 

The following restrictions on the xt-y 
and xl^ parameters were imposed: xl-y < 
25, X37 > 20. The first restriction caus- 
ing only 5% loss of actual rjn^j 
events reduces background from the pro- 
cess by approximately 30% and al- 
most completely removes background from 
the process (||). The second cut sup- 
presses background from the processes 
— ^ 777 37, (p — > 7r°7 — > 37, e+e^ — > 
27,37 (QED). 

For further background suppression, 
an event configuration (the photon en- 
ergies and angles after 5-7 kinematic fit- 
ting) was compared using modification 
of kinematic fitting technique developed 
in the work with ones expected for 
the process (||) and background processes 
@, (H). Corresponding measure of dis- 
crepancy P is an increase of for a 
5-7 event after application of additional 
requirements on intermediate states for 
each tested hypothesis. The following 
hypotheses were considered: 



• an event is a cascade reaction e^e^ - 
X7, X rjTT^ where X is some 
intermediate particle; the P^^^ pa- 
rameter and invariant masses of 
photon pairs, presumably produced 
in the decays of vr*^ and rj mesons 

( M,r and ) were calculated; 

• an event is a cascade reaction e^e^ - 
X7, X 7r°7r°; the P^^r^ param- 
eter was calculated; 

• an event is of the process e~^e~ — 
ci;7r°,u; 7r°7; parameter M^^ — 
an invariant mass of n^'-f pair from 
the u TT^'-f decay, was calcu- 
lated; 

Relative contributions from the back- 
ground processes (H) and (H) vary with 
m^TT — invariant mass of r^vr*^ pair. At 
rrirjTT < 975 MeV the background from 
the process @ becomes significant. Ad- 
ditional cut PTrn-y > 2 supprcsscs it by a 
factor of three reducing detection effi- 
ciency for the process (0) by only 12%. 
At rrirjTT < 900 MeV the dominant back- 
ground comes from the process (|^). In 
this region restriction < 725 MeV 
removes background almost completely. 

The scatter plot for invariant mass 
versus max / Ebeam — normalized en- 



ergy of the most energetic photon in the 
selected events is shown in fig.|l], where 
two regions are distinguishable: 

max / Eiyeam > 0.68 dominated by back- 
ground from the reaction with a nearly 
uniform distribution and 
E^max/Ebeam < 0.68 whcrc the back- 
ground is small and the points are grouped 
close to 77-meson mass. The vs. M^^o 
distribution in events with E^max/ Ebeam < 
0.68 is shown in fig.^. It is clearly peaked 
at 7r° and 77-meson masses. 

For final selection of — > rj'K'-f events, 
in addition to the cuts described above 
the restriction E-ymax/ Ebeam < 0.68 prac- 
tically completely removing background 
from the process (^) and PrjTTj < 7 were 
applied. The Pr)7r7 distribution for the 
selected events is shown in fig.^ The 

CU.tj -Pfjj^-y 7 roughly corresponds to 

the restrictions | M^o — m^o | < 30 Me V 
and — < 30 MeV, where m^o 
and rrirf are the 7r° and 77 masses. Total 
of 39 events were found with expected 
background of 3.2 ± 0.7 events. In the 
region 7 < Pj^t,^ < 15 ten events were 
found in agreement with estimated 8.9± 
0.4 ?77r°7 events plus 6 ± 1 background 
events. Thus, after cuts described above 
the event sample still contains background 



of about 10%. After background sub- 
traction 35.8 ± 6.3 events of the process 
e~^e~ riTT^'-f are left. 

Data analysis. Fig.| shows cos a 
distribution for the selected events, where 
a is an angle between the recoil pho- 
ton in the reaction (|^) and r^-meson mo- 
mentum in the rjir^ rest frame. Esti- 
mated background of 3.2 events is sub- 
tracted. Experimental distribution is in 
a good agreement with the simulated 
one, which was initially isotropic as ex- 
pected for a scalar intermediate state. 
Its visible slope is a consequence of the 
E-ymax < 0.68 cut. Such an agreement 
may be considered as an evidence that 
?77r°-system is produced in a scalar state 
{P{x'^) = 61%). In fig.| the cos^^ dis- 
tribution is shown. The 6^ is a polar an- 
gle of the recoil photon in the reaction 
(g). It also agrees (P(x^) = 22%) with 
the simulated distribution (1 + cos^ 6^/) 
expected for production of a scalar par- 
ticle and photon. 

Detection efficiency obtained by sim- 
ulation must be corrected for event loss 
due to additional spurious photons and 
for imprecise simulation of parameters 
used in the event selection cuts. Corre- 
sponding correction factor was obtained 



from experimental data. To this end 
cross section of the process (H) was mea- 
sured using the selection criteria similar 
to those described above for the process 
@. The result was compared with our 
earlier measurement [1^. It was found 
that simulation overestimates detection 
efficiency for the process (0) by 5%. 

In the Table |I] the numbers of se- 
lected events, detection efficiencies, and 
measured differential branching ratios dB/dm 
as a function of m^^r invariant mass are 
listed. The detection efficiencies and 
dB / dm values are given at middle points 
of the corresponding invariant mass bins. 
Uniformly distributed background of 3 
events was subtracted. The detection 
efficiency averaged over the experimen- 
tal invariant mass spectrum is equal to 
2.1%. 

Fitting of energy dependence of the 
experimental cross section by the sum of 
the resonant cross section of the process 
(H) and energy-independent background 
results in 

5(0 ^ r/7r°7) = (0.88±0.14±0.09)■10-^ 

(7) 

The main sources of systematic error 
here are uncertainties in the measured 



cross section of the process and in 
average detection efficiency due to large 
statisical error of the observed vjii^ in- 
variant mass spectrum. 

In fig.|] the dependence of the mea- 
sured — > r77r°7 decay branching ratio 
on the invariant mass of the r^vr^ pair is 
shown. In spite of smaller recoil pho- 
ton phase space at high rjn^ invariant 
masses the observed mass spectrum shows 
enhancement in this region. This means 
that r/vr^ system is produced in some 
resonant state. The only known res- 
onance which have relevant mass and 
quantum numbers is ao(980) and the ob- 
served enhancement at large ?77r° invari- 
ant masses can be described as manifes- 
tation of ^ ao7 decay. Known — >■ 
p^7T^,p —>■ rj'y decay mechanism must 
produce rjTc^ pairs with smaller invariant 
masses and as was already mentioned 
its branching ratio is much smaller than 
observed one, although its amplitude should 
be taken into account in the approxi- 
mation of the whole mass spectrum in 
future high statistics experiments. For 
M^^ > 900 MeV we have: 

5(0 ^ 777r°7) = (0.46±0.13)-10"^ (8) 



Discussion. Since the experimen- 
tal data show large contribution from 

the (f) — > ao7 decay an attempt was made 
to approximate observed invariant mass 
spectrum in assumption of pure (p 
ao7 and assuming decay dynamics as 
described in the work 0. This hypoth- 
esis gives rather good approximation of 
the experimental data. The fitting curve 
is shown in fig.| (-P(x^) = 61%). The 
following optimal values of the Oo-meson 
parameters were obtained: 

Ma, = 995^f^MeV 

9lK^K-/^^ = {^-^^-lt)GeY' (9) 

9l,J^^ = (0.77tS;i)GeV^ 

The corresponding fitting curve is shown 
in fig.p. The optimum value of the ratio 
9aov7r/9aoK+K- = 0.75t[J;i wlthlu exper- 
imental errors satisfies the relation be- 
tween coupling constants gaorin = 
9aoK+K- obtained in 0| in the assump- 
tion of 4-quark structure of meson. 
If we fix this ratio at its 4-quark model 
prediction, the optimal values of other 
fit parameters become: 

Ma, = 994lf MeV ^^^^ 
^^„K+x-/4^ = (l-05li)GeV^ 



The mass Ma, is in agreement with the 
PDG value 983.4 MeV fl^. 

Comparison with (H) shows that about 
50% of the observed branching ratio (0) 
corresponds to M^^ > 900 MeV and 
the observed invariant mass spectrum is 
consistent with the model 0. Thus it 
can be assumed, that the ao7 intermedi- 
ate state dominates in this decay. Other 
decay mechanisms may contribute, for 
example P^t^^,P Vli although 
rough estimation shows, that its contri- 
bution can be neglected at present level 
of experimental errors (|^). Assuming 
pure — > ao7 decay, we get 

5(0 ^ ao7) = (0.88±0.17)-10"^ (11) 

Conclusions. In this work, us- 
ing experimental data corresponding to 
about 2 ■ 10'' produced mesons, 36 
events of the rjir^'y decay were found. 
The measured branching ratio of this 
decay fi(0 ^ r^vr^) = (0.88 ± 0.14 ± 
0.09) • 10~^ is in agreement with our pre- 
vious result -8(0 ?77r°7) = (0.83 ± 
0.23) ■ 10^^ 0, based on analysis of a 
part of the experimental statistics, as 
well as with the CMD-2 measurement: 
5(0 ^ r/7r°7) = (0.90±0.24±0.10)-10-^ 



0]. Observed enhancement in the rjir^- 
pair invariant mass spectrum at large 
masses shows large contribution of ao7 
intermediate state. Assuming dominance 
of this mechanism we obtain B{(f) —>■ 
ao7) = (0.88 ±0.17) ■ 10"^ 
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Figure captions 

• Figure |I]: Distribution of — 
reconstructed mass of rj meson ver- 
sus the energy of the most ener- 
getic photon in the event E^rnax/ Ebeam- 

• Figure 0: Distribution of — 
reconstructed mass of rj meson ver- 
sus — reconstructed mass of 

7r° for events with E^rnax/ Ei,eam < 

0.68. 

• F2(7nre ^ The P^^^ distribution. Points 
with error bars - experimental data. 
Histogram — simulated signal from 

— > ?77r'^7 decay corresponding to 
branching ratio of 0.9-10"'^, shaded 
histogram — estimated background 
from the e'^e^ ujtt^ and — >■ 
777,/o(980)7 processes. 

• Fi^rwre^ The cos a distribution, a 
is an angle between recoil photon 
and rj meson in the r/vr" rest frame 
for selected r]TT^'^ events. Points 
with error bars - experimental data, 
histogram - simulation of the pro- 
cess (|T]) with BR = 0.9 ■ 10^^ 

• Figure ^ Recoil photon polar an- 
gle distribution for selected rjn^'j 
events. Points with error bars - 



experimental data, histogram - sim- 
ulation of the process with BR = 
0.9- 10-^ 

Figure^ Mjjn invariant mass spec- 
trum. The fitted curve corresponds 
toM,„ = 995tll MeVgl^+j,./4n = 
{lAtli)GeV' gU^n = (0.77lJ:l)GeV^ 




0.6 0.7 0.8 0.9 

^y.ma:/^ beam 



Figure 1: Distribution of — re- 
constructed mass of 7] meson versus the 
energy of the most energetic photon in 
the event -E^^m^./ Ei^g^Q^jy^. 




Figure 2: Distribution of — recon- 
structed mass of 7] meson versus M^r — 
reconstructed mass of 7r° for events with 

E-y max / Efyga^ < 0.68. 




Figure 3: The P^^^ distribution. 
Points with error bars - experimental 
data. Histogram — simulated signal 
from ?7'7r°7 decay corresponding to 
branching ratio of 0.9 ■ 10~^, shaded his- 
togram — estimated background from 



the e~^e 
processes. 



cu7r° and (f) r/7,/0 (980)7 




Figure 5: Recoil photon polar an- 
gle distribution for selected rjn^'y events. 
Points with error bars - experimental 
data, histogram - simulation of the pro- 
cess (I) with BR = 0.9 ■ 10"^ 
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Figure 4: The cosa distribution, a is 
an angle between recoil photon and 77 
meson in the r]7T^ rest frame for selected 
r77r°7 events. Points with error bars - ex- 
perimental data, histogram - simulation 
of the process (|l|) with BR = 0.9 ■ 10"^. 
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Figure 6: 
spectrum. 



MriTT invariant mass 
The fitted curve cor- 
responds to Mao = 995lio MeV 
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Tabic 1: The number N of found (j) 
r77r°7 decay events, detection efficiency 
e, and measured dB/dm as a function 
of rrinn — the tjti^ invariant mass. 
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